Fidelity of qubit rotations in the presence of thermal motion 
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We investigate both theoretically and experimentally the effect of real-life imperfections on the 
fidelity of qubit rotations in the presence of thermal motion. Specifically we develop a model 
predicting the coherence time of Rabi oscillations in the presence of an inhomogeneous driving 
field or Ramsey oscillations in the presence of an inhomogeneous state-selective potential. We 
demonstrate the theoretical dependence by an experiment observing the coherence time in a volume 
much smaller than the entire sample. 
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Since the middle of the previous century, the two-level 
system has been the source of numerous new theoretical 
and experimental insights into the basic tenets of quan- 
tum theory. Later on, the two-level system has become 
the main building block of numerous applications rang- 
ing from clocks [1 and magnetic sensors [2] to quantum 
communications 3 and computing [JJ. Recently, this 
two-level system has been named the qubit. Understand- 
ing the sources of errors when manipulating this system is 
of paramount importance for technological applications 
as well as fundamental studies [5]. 
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FIG. 1: (color online) (a) absorption image of S7 Rb atoms 
at the F = 2 energy level after applying a spatially inhomo- 
geneous MW field (T = 43 /LtK). The image is taken while 
in free fall and the observed spatial axes are perpendicular 
to the gravitational axis. The MW horn antenna (depicted) 
radiates electromagnetic waves with decreasing amplitude as 
a function of distance, thus generating along the cloud mul- 
tiple Rabi oscillations which can be viewed simultaneously, 
thereby creating a fringe-like pattern, (b) Rabi oscillations 
measured 10 cm away from the antenna [e.g. along the ver- 
tical line in (a)]. Each data point represents the same pixels 
in the cloud at different MW pulse durations. The coherence 
decay is due to the gradient of the MW field and the velocity 
of the atoms (temperature), (c) The plot shows the depen- 
dence of the observed Rabi frequency on the position across 
the cloud where the origin of the graph refers to the center of 
the cloud. The data was fitted to a linear curve and results 
with d x Q, = -1.742 ±0.129 [1/msec ■ mm]. 



The manipulation and evolution of such a two-level 
system, defined in quantum theory as \ip) = a\l) + b\2), 
is perhaps best described by a vector on the Bloch 
sphere defined by angles and <f>, where a — cos | and 
b = sin|e 11 ^. In a stationary situation the two qubit 
states may be chosen as energy eigenstates and then 
4> = —{E% — E\)t/h (Ei being the energy of state i). 
Any manipulation of the qubit state is typically called a 
rotation as it rotates the state vector on the sphere. For 
example, in the presence of an oscillating driving field 
resonant with the energy E2 — E\ the angle 6 rotates 
as 6 = fit, where the Rabi frequency Q is proportional 
to the amplitude of the driving field. The fidelity of a 
single-qubit operation is defined as the overlap of the ac- 
tual final state with the target state. 

In real-life devices, gradients of static fields (determin- 
ing E 2 — Ei) or driving fields unfortunately exist (field 
gradients are an experimental reality even if cavities are 
used [HE])- In such a case, there exists a first order 
broadening effect which reduces the coherence time. This 
may be viewed in Fig. [l|a), where a gradient of a driving 
micro- wave (MW) field induces multiple Rabi frequencies 
in space. This is true also for a sample temperature of 
T = 0. However, when T ^ a second order deco- 
herence effect appears, as atoms move between different 
Rabi frequency regions. This decrease may be observed 
by studying the Rabi oscillations in a small volume cle- 
ment, much smaller than the sample size [Fig. [ljb)]. It 
is this second order effect which is the topic of this paper. 

We investigate both theoretically and experimentally 
the reduction in fidelity in a small volume element as a 
function of atom temperature and field gradient. The 
significance of the small volume element fidelity goes be- 
yond fundamental interest; it may also be of practical 
importance, for example, in the optimization of signal to 
noise in situations in which the coherence time is more 
dominant than the number of atoms, or when, for techni- 
cal reasons, imaging a volume smaller than the trap size, 
e.g. when utilizing near field optics as enabled by atom 
chips [81-11 1 j . Next, we calculate the fidelity of the com- 
mon rotation of 6 = 90°, typically called a ir/2 pulse. 
Finally, we show that our model may be expanded to 
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describe any trapped system in which both qubit states 
feel the same potential (e.g. \F,m,F) = |1, —1) and |2, 1) 
[14] ) . or systems with arbitrary qubit states where a con- 
stant potential is used i.e. box potentials (e.g. [T5]). 
Such traps may be used in the future for a wide variety 
of applications (e.g. quantum computers) as they reduce 
potential induced inhomogeneities. Our results may be 
used to improve the fidelity of such systems. 

The model system we are analyzing in this paper is 
that of freely propagating atoms when no potential is 
applied. We first construct a theoretical model which 
is solved both analytically and numerically, and then ex- 
perimentally demonstrate the consequences of the model. 
Our experiment utilizes a far field MW imaging method 
presented here for the first time [T^] [Fig. [TJa)]. A sim- 
ilar near-field method was developed independently by 
the Treutlein group [13] , 

Our theoretical model assumes an ensemble of two- 
level atoms whose motional degrees of freedom evolve 
classically in time and are given initially by a classical 
distribution P(x, v) in position and velocity space. The 
internal state of an atom with an initial position Xo and 
velocity vo is generally given by 

IK,v„(*)> = Oxo,v (*)|l) + e- iwi2t & Xo ,v (t)|2), (1) 

where uj\ 2 = (E2 — E\) /h is the transition frequency. The 
equations of motion for the evolution of the internal state 
under the influence of an applied radiation field are 

«x„,v = -^-" 12 ^[x(x ,v ,i)]6 Xo , Vo 

ix„,v = -ie-^-^^nixfxo.vo,*)]*^^, (2) 
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FIG. 2: (color online) Rabi oscillations for different temper- 
atures (8 fjK to 37 /iK). The graphs present the normalized 
population in F — 2 as a function of the MW pulse duration. 
It can be seen that the coherence time of the oscillations in- 
creases when the temperature decreases. Each of the graphs 
was fitted to the model of Eq. ( 10 1 (see text for more details). 



where uj is the frequency of the driving field, and Vl (x) 
is the Rabi frequency for qubit rotation at a given point 
x(xo,vo,t) along the atomic trajectory, which is propor- 
tional to the amplitude of the driving field at that point. 
As the driving field is usually in the microwave frequency 
regime, we may ignore the Doppler shift and assume a 
constant frequency for all velocities. 

Once Eqs. Q are solved for every trajectory, the den- 
sity matrix for the internal state is given by an integral 
over the initial external degrees of freedom of the atoms 

Pij(x,t) = J rf 3 x y d 3 v P(x ,vo)(5 3 [x- x(x ,v Q ,t)] x 



X Wxo.voWXV'xo.voWlj) 



(3) 



for i,j = 1,2, such that (1|V>) = a and (2\tfj) = b. 

A simple analytic solution for Eqs. ^ may be obtained 
for the case where the driving field is resonant with the 
atomic frequency uji 2 . Then the Rabi frequency fi may 
be chosen to be real. In this case we obtain the solution 



a(t) = cos[£(t)]a(0)-isin[£(i)]&(0) 
b(t) = -ism[£(i)]a(0) + cos[£(t)]6(G), 



(4) 



where £(t) = § fj. \dt'VL(t'). 

First we consider the simple case which applies to 
the conditions of our experiment, as described below. 
We assume a free collisionless motion during the time 
where the radiation is applied, such that the velocity v 
of a particle remains equal to the initial velocity Vg and 
x(xq,v, t) = xo + vt. We also assume that the distances 
|v|t traversed by the atoms are smaller than the size of 
the initial atomic cloud. We further assume a small field 
variation such that over these distances 



n(x) ft + (vn) • (x - x c ) 



(5) 



where Vi7 is a constant gradient and fio is the Rabi fre- 
quency at a given central point x c . In this case the phase 
£(t) in Eq. @ is 

e v (x,f) = ^2(x)t- Jv-(Vr!)i 2 , (6) 

where fl (x) is the Rabi frequency at the observation 
point. As only the velocity component along the di- 
rection of the field gradient affects the internal atomic 
dynamics, we take this direction to be denoted by x and 
refer to v x as v. Assuming a Boltzmann distribution of 
the atomic velocity 

™„,2 



P(v) 
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exp 



2k B T 



(7) 



The resulting t v and \ are shown. 



and taking the atoms to be initially at state |1), we find, 
using Eq. ([3]), that the population of the state |2) is 
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dv P(w)sin 2 [^(a;,t)] 
[1 - cosO(a;)texp[-t 4 (9 x r2) 2 fc B T/8m] 
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where we have used sin 2 £ = (1 — cos2£)/2. It fol- 
lows that the amplitude of Rabi oscillations decays as 
exp[— (t/r v ) 4 }, where the temperature dependent coher- 
ence time is 



(v\d x n\y/*' 



(9) 



with v — \J (v 2 ) = y^ksT/m being the average velocity. 

If the imaging system which probes the state of the 
atoms at a given point has a finite resolution, such that 
every pixel of the camera represents a convolution of the 
atomic density at a given point with a Gaussian function 
(V^ cr i)~ 1 exp[—(x~x') 2 /2aj] of radius 07 due to imaging 
limitations, then the number of atoms at level \2) is 



n 2 (x,t) = - {1 -cos[n(x)t]exp[-{t/T x ) 2 - (t/r v ) 4 }} 



where 



2 i/2 



(10) 

(11) 



This integration over a finite spatial range represents the 
ordinary decoherence effect of the spatial inhomogeneity. 

To demonstrate the theoretical model we experimen- 
tally investigate a cloud of freely propagating atoms (in 
free-fall). We start our experiment with a cloud of 10 6 
cold 87 Rb atoms. The atoms are cooled by a standard 
Magneto-Optical- Trap (MOT) followed by laser molasses 
to the required temperature, down to a few /zK. The 
atoms are prepared in the F = 1 state. To avoid in- 
homogeneities, the atoms are then released into free-fall 
for the duration of the experiment, typically 10 — 30 ms. 
During the free-fall the atoms are subjected for a time t, 
via a horn antenna, to a MW field tuned at 6.8 GHz to 
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FIG. 3: (color online) A comparison of fitting different mod- 
els to the data (T = i3(iK). The dotted, dashed and solid 
lines present the envelope of the exponential (t) mo del, the 
Gaussian (t 2 ) model and our model presented in Eq. ( 10 1, re- 
spectively. The latter model returns a x 2 value of 0.972, while 
the Gaussian model returns \ 2 ~ 0.947 and the exponential 
model returns \ 2 = 0.844. 



the \F,m F ) = ]1, 0> = |1) |2,0) = |2) transition. After 
time t the MW field is switched off and the population 
of the atoms in the F = 2 state is determined by on 
resonance absorption imaging directed along the gravita- 
tional axis. The MW field induces Rabi oscillations, the 
coherence time of which we analyze. 

In Fig. [lja) we present a typical experimental image 
following the sequence described above. As the horn an- 
tenna produces a spatially inhomogeneous MW field, a 
gradient of Rabi frequencies is produced along the cloud 
(decreasing frequency with growing distance from the an- 
tenna) ; this is exhibited in the figure as a fringe- like pat- 
tern [12"] . In Fig. [ijb) we present a typical Rabi os- 
cillation plot (F = 2 population vs. t) as seen in the 
data taken from a vertical strip of pixels one pixel wide 
and 100 pixels long (perpendicular to the direction of the 
field gradient). One may deduce the gradient by fitting 
the spatial dependence of the measured Rabi frequencies 
to a linear slope. Far away from the antenna one may 
assume a nearly linear slope as one expects the magnetic 
field radiated by the antenna to fall as 1/r [Fig. [TJc)]. 

Next, we analyze the coherence time as a function of 
temperature. In Fig. [2] we present as an example four 
data sets of Rabi oscillations at different temperatures. 
The model presented in Eq. (10 1 fits well to the data. 
To make the fit, we insert into Eq. (10) an arbitrary 



constant phase to account for possible systematic shifts 
in the timing of the experiment, so that the fit is to the 
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FIG. 4: (color online) Coherence time as a function of the 
temperature (T = 8 — 102 ^K). The data points are the co- 
herence time extracted from the t 4 model, fitted to the data 
as shown in Fig. [2] The temperatures were measured using 
time-of flight (TOF). The error values of the coherence time 
and temperature were calculated in the model with a confi- 
dence level of 95%. The data points of coherence time (t„) for 
each temperature (T) were then fitted to the model presented 
in Eq. S and the fit returns a x 2 value of 0.996. The region 
between the dotted lines indicates the range of the predicted 
theoretical values taking into account the errors on the mea- 
sured Rabi frequency gradient. This confirms to a high level 
of confidence our theoretical model. 
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function /(t) = Aexp[-(t/r x ) 2 -(t/r v ) 4 ] cos(fit+0)+£. 
The A and i? constants are added because the oscillation 
amplitude and the average are not exactly known. In 
addition, to make the fit, we first need to estimate t x . 
When leaving both t x and t v as free parameters, a fit on 
the T — 43 fiK data set appearing in Fig. 1(b), returns 
t x = 8.8 msec with a \ 2 of 0.97. This corresponds to 07 of 
94 microns [Eq. (11 1, with d x fl = 1.7 (mm msec) -1 from 
Fig. 1(c)] which in turn corresponds to a misalignment 
of our 30 cm focal length lens by 1 mm or so along the 
imaging axis. As we estimate that our optics alignment 
error is at least that (as the cloud size itself is about 

1 mm in all directions), we adopt in the following the 
t x = 8.8 msec value [Eq. ( 11 )] in the rest of the paper and 
leave only t v as a free parameter. Let us note that the 
procedure is very robust and changing 07 by a factor of 

2 in each direction returns \ 2 values with a mere change 
of 1%. Finally, we also fit the data to Eq. (10) while 



replacing the power of 4 by a free parameter d and find 
that it converges to values of d — 3.778 — 4.1 with \ 2 
values 0.961 - 0.971. 

For comparison we plot in Fig. [3] one set of Rabi oscil- 
lations (T = 43 /iK) with a fit to three possible models: a 
Gaussian model, an exponential model and the t 4 model. 
One can clearly see that the t 4 model better fits the data 
by far. We now use the gradient measured by the fit pre- 
sented in Fig. |ljc) to compare the observed coherence 
times at different temperatures to the theoretically ex- 
pected value. As presented in Fig. |4j the data fits well 
with the theoretical prediction. 

The t 4 exponential dependence of the decay of Rabi 
oscillations follows from the Gaussian velocity distribu- 
tion in the thermal cloud. Atoms with higher velocity 
travel a larger distance along the field gradient thereby 
acquiring a larger phase difference relative to atoms at 
rest in the detection point. This additional phase is an 
integral over time along the way which was traversed by 
an atom with a given velocity, which is vt, such that the 
total phase depends quadratically on time. The combi- 
nation of the quadratic dependence of the phase on time 
and the quadratic exponential dependence of the distri- 
bution on velocity provides the t 4 quadratic dependence 
of the coherence on time. We may consider t v as a criti- 
cal time such that at smaller times the qubit rotation is 
unaffected by the velocity. 

Finally, we use the t 4 model to infer the fidelity of a 
7r/2 pulse. Fidelity is defined by the overlap between 
a target state ^target and an actual state \ip). If the 
actual state is not pure then it is described by a den- 
sity matrix p = Wj\il)j){il)j\. The fidelity is then 

given by F = £\ w 3 \(ipA tp) targct | 2 . For a ir/2 pulse 

the target state is |i/7targct = cos(7r/4)|l) — i sin(7r/4)|2). 
For a given velocity the actual state is \ip v ) — 
cos (|(7r + d x CLvtQ)] |1) — ism [|(tt + dxClvt 2 )] |2), where 
to = 7t/457q. The overlap between the actual state 



and the target state is {ip v \tp) target = cos (?; 9a; 0^/4). 
Integrating the square of the overlap over the differ- 
ent velocities we obtain F 2 = J dvP(v)\(ip v \i/jtar g et)\ 2 = 
\ {1 + exp[-(7r/4fi r t ,) 4 ]}. When Q t v > tt/4 the fi- 
delity is almost 1, while if Q,ot v < tt/4 the fidelity drops 
to a minimum value of F = l/y/2, which represents the 
fidelity for a totally random actual qubit state. 

Note that an analysis similar to the one presented 
above for the effect of a gradient of a driving field applies 
also to a gradient of a stationary field when no driving 
field is applied. Such a stationary field consists a local 
differential potential V(x) for the two levels, such that 
their mutual phase cj) has a local term </>(x) = — V(x)t/H. 
In this case the local coherence between the two levels is 



p 12 (x) = a(0)&*(0)e 



d 3 vP(x,v) x 



x exp 



diV[x(0] 



(12) 



resulting in coherence times similar to t v and t x above, 
with the replacement Vf2 — > VV/h. This result makes 
the analysis above relevant to situations such as the dark 
time in a Ramsey process in, for example, atomic clocks. 

As an outlook we note that our model is relevant also 
to systems in which both qubit states feel the same po- 
tential, or in which for arbitrary qubit states the poten- 
tial is constant. As is many times the case, when the 
position dependent interaction with the environment is 
stronger than the energy dependent interaction |16j and 
consequently the position base coherence time is smaller 
than the energy base coherence time, potentials are ex- 
pected to house classical particles in the sense of external 
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FIG. 5: (color online) Numerical simulation: Coherence of 
Rabi oscillations in the middle of a box trap with length L 
and a driving field gradient along the same axis, as a function 
of velocity spread Av [units of L/(arb. time units)]. The 
gradient is d x fl = 0.1 [units of 1/(L- arb. time units)]. 
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degrees of freedom and quantum particles as far as the 
internal degrees of freedom are concerned |17j . Thus our 
model which is based on classical dynamics should also 
apply to trapped particles (as well a single trapped par- 
ticle). As an example we analyze in the following a box 
potential (for which the above model for freely propagat- 
ing atoms may be taken as the limit of a very large box) . 
To adapt our model to box potentials one has to take into 
account the reflections at the walls. Our numerical study 
which will be presented elsewhere shows that indeed the 
coherence initially drops like exp (— av 2 t A ) at short times 
such that t < L/v (where L is the size of the box). For 
longer times, the addition of the walls adds the following 
behavior: The coherence time reaches a minimum value 
and oscillates around this value with decreasing oscilla- 
tion. Eventually it stabilizes around a constant value 
which grows with velocity and is high when \v\t v 3> L 
(t v being the coherence time for L — > oo). This is pre- 
sented in Fig. [5] Higher velocities give higher coherence 
for longer times because higher velocities correspond to 
paths with multiple reflections so that the field gradient 
is averaged out. This may be understood by making an 
analogy between the reflections at the walls and multi- 
ple 7r echo pulses in the middle of a Ramsey sequence 
[18]. An analogy can also be made to motional and col- 
lisional narrowing |19| . Hence the walls in fact simply 
add motional or collisional narrowing. However, if the 
population of a certain qubit level over the whole area of 
the box is probed, then the coherence functions in Fig. 
5 must be multiplied by the spatial integration factor as 
in Eq. (10), which is, in our case, s'm(d x ilLt)/(d x flLt), 



corresponding to a decoherence due to ordinary inhomo- 
geneous broadening t x = (dxilL)^ 1 . 

To conclude, we have shown that gradients of DC or 
AC fields reduce the coherence of qubit rotations both di- 
rectly though an inhomogeneous rotation frequency and 
through a second order temperature dependent effect due 
to particle motion. We have demonstrated our theo- 
retical model by a simple experimental system of freely 
falling atoms and also by a numerical example of a par- 
ticle in a box. Our results may be useful to improve the 
operation of quantum devices. 

We gratefully acknowledge help from the members of 
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mon Machluf, Menachem Givon and Zina Binshtok. 
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